Ab initio GW electron-electron interaction effects in Quantum Transport 
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We present an ab initio approach to electronic transport in nanoscale systems which includes 
electronic correlations through the GW approximation. With respect to Landauer approaches based 
on density-functional theory (DFT) , we introduce a physical quasiparticle electronic-structure into a 
non-equilibrium Green's function theory framework. We use an equilibrium non-selfconsistent G^VF" 
self-energy considering both full non-hermiticity and dynamical effects. The method is applied to 
a real system, a gold mono-atomic chain. With respect to DFT results, the conductance profile is 
modified and reduced by to the introduction of diffusion and loss-of-coherence effects. The linear 
response conductance characteristic appear to be in agreement with experimental results. 
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Electronics at the nanoscale, namely nanoelectronics, 
represents the next years' technological challenge. It 
is boosted not only by the need for shorter integration 
scales, but also by the expectation that unusual quan- 
tum effects [H are going to be observed due to quan- 
tum phenomena effects. Beside the experimental efforts 
to synthesize nanoelectronic devices, quantum transport 
theory Q has the formidable task to understand and to 
model the mechanisms behind these phenomena and to 
predict them from a first principles approach. 

In the last years, a combination of ab initio density- 
fimctional theory (DFT) calculations together with 
the description of transport properties in a Landauer- 
Biittiker (LB) framework j3| has demonstrated its ability 
to describe small bias coherent transport in nano junc- 
tions 0, 0j 0| • These approaches were successful in ac- 
counting for the contact resistance and conductance de- 
grading mechanisms induced by impurities, defects and 
non-commensurability patterns in the conductor region. 
The major objections raised to such method are: (z) 
the Kohn-Sham (KS) electronic structure is in principle 
unphysical, to be considered only as an approximation 
to the quasiparticle (QP) electronic structure; (ii) non- 
coherent and dissipative effects due to electron-phonon 
{e-ph) and electron-electron (e-e) scattering can be taken 
into account only approximatively in the LB formalism; 
(Hi) non-linear response and far from equilibrium finite- 
bias transport are not accessible, since DFT cannot be 
applied to open systems and is not a non-equilibrium the- 
ory (although recent works 6] have demonstrated that 
time- dependent DFT can tackle the problem). 

Non-equilibrium Green's function (NEGF) theory [7, 
[§] is in principle a correct approach to address the above 
objections. The critical point within this theory is the 
choice of good approximations to the self-energy E*", and 
coherently to the scattering functions S^'^. This en- 
sures that both the renormalization of the QP energies 



and the electron diffusion mechanisms due e.g. to e-ph or 
e-e interactions will be properly taken into account. First 
works studying the role of the e-ph coupHng and of 

short-range e-e correlations [ill , [l^ or the renormaliza- 
tion of the QP energies [l3| have recently appeared in the 
literature, or, we are aware, are going to appear [l^ . The 
role of correlations, apart from being central in explain- 
ing e.g. Coulomb blockade and Kondo effects, could also 
be crucial in bridging the gap between experimentally- 
observed and LB-predicted conductances, in some cases 
are orders of mag nitude off @,[li,[il. 



In this work, we introduce electronic correlations in 
the calculation of transport by an ab initio approach 



based on Hedin's GW approximation (GWA) [1 
in the framework of NEGF. In our scheme, the GW 
self-energy is built at equilibrium and the Green's func- 
tion is calculated by direct solution of the Dyson equa- 
tion. For the lead/conductor/lead geometry, the GW 
self-energy is summed to the lead's self-energies; the 
electronic conductance is calculated through the Meir- 
Wingreen formula [3], a NEGF Landauer-like expres- 
sion derived for interacting conductors. We apply this 
scheme to a realistic system, a gold mono-atomic unidi- 
mensional chain [i^] , and we study the effects induced on 
transport properties by the different components of the 
GW self-energy, the hermitean and the non-hermitean 
parts and the dynamical dependence. Our results show 
that the conductance profile is considerably modified by 
the real part of the GW correction. The imaginary part 
introduces a suppression of the conductance, which is 
negligible close to the Fermi energy, but that increases 
with the energy. Finally, the full dynamical dependence 
of the GW self-energy introduces further structures far 
from the Fermi energy, which have to be ascribed to satel- 
lite excitations of the system. The GW smooth drop on 
the conductance characteristics as a function of the bias 
at very low voltage compares favorably with the trend 
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FIG. 1: 2-particle Green's function in the self-consistent 
Hartree-Fock (a) and in the non self-consistent G'^W^ ap- 
proximation (b). Thin line: non self-consistent G"; double 
line: self-consistent G; wiggly line: RPA non self-consistent 



experimentally measured in gold nanowires [20|. 

With respect to Hartree-Fock (HF), which already 
renormalizes the energies for the e-e classical repulsion 
and exchange, the GWA introduces the important con- 
tribution due to correlations and to e-e scattering dif- 
fusion mechanisms responsible for loss-of-coherence in 
transport. Indeed, by direct inspection of the diagram- 
matic representation of the 2-particle Green's function 
G2 (see Fig. [T] and Refs. 0, [ll]), one can see that the 
G^^ describes an uncorrelated propagation and that col- 
lisional terms are missing. This implies that the HF scat- 
tering functions, and S^, are exactly zero. On the 
other hand, even a non-selfconsistent G'^W'^ approxima- 
tion introduces a coUisional term [the last diagram in 
Fig. [Tfb)] which gives rise to non-zero S^''', and in turn 
to e-e scattering mechanisms and incoherent, dissipative 
effects in transport. As it is shown by the corresponding 
G2 Feynman diagram, the G°M^° approximation is not 
a conserving approximation in the Baym and Kadanoff 
sense 0, leading to, e.g., non conservation of the num- 
ber of particles. However, the relative deviation from 
the exact density brought by the G°PF° approximation 
has been evaluated by Schindlmayr et al. [18f| to be only 
of the order of 0.05% for the range of metallic densities 
(r^" = 3.01) of interest here. 

Our starting point is a standard DFT-LDA calcula- 
tion based on plane waves (PWs) and norm-conserving 
pseudopotentials for an infinite mono-atomic chain of 
gold atoms using periodic boundary conditions 2]J . The 
KS electronic structure is calculated both at the relaxed 
atomic distance (4.72 Bohr) and in a stretched geome- 
try (5.32 Bohr), so as to simulate the experimental sit- 
uation described in Ref. ^2fl,] (conductance measures of 
a gold monoatomic chain pulled up from a gold surface 
by an STM tip) and also the calculations reported in 
Ref. 0. ^From the DFT KS eigenfunctions, we obtain 
an orthonormal set of maximally localized Wannier func- 
tions (MLWF) [22I, which are used as a basis set in the 
calculation of quantum transport. The following step is 
a converged [23| GW plane- wave calculation of both the 
QP energies and the self-energy matrix elements for the 
six bands around the Fermi energy, corresponding to the 



gold sd-manifold. The self-energy in the G^W^ approxi- 
mation at equilibrium is given by: 

/oo 
dLo'e-'^'^^G\i^-Lo')W\Lo'), (1) 
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where G° is the Green's function built on the non- 
interacting KS electronic structure and is the dy- 
namically screened interaction given by the RPA polar- 
izability, p^^^ = —iG^G^. Since for transport we need 
a fine-grid fully dynamical dependence of the self-energy, 
we calculate the frequency integral of Eq. ^ in three dif- 
ferent ways: (z) by approximating the dynamical depen- 
dence of W{uj') through a plasmon-pole (PP) model 17|; 
(a) by a contour deformation (CD) method [2^ . which 
consists in a deformation of the real axis contour such 
that the self-energy can be calculated as an integral along 
the imaginary axis minus a contribution arising from the 
residual of the contour-included poles of G; (Hi) by an 
analytic continuation (AC) method [25] . i.e. calculating 
the integral and also the self-energy on the imaginary 
axis and then performing an analytical continuation to 
the real axis. In the last step we carried out the quan- 
tum transport calculation using a modified version of the 
Want code d El- We first projected the GW self- 
energy, as well as the non- interacting hamiltonian H^, 
on the Wannier functions basis set (non-diagonal self- 
energy elements in the Bloch basis were neglected). We 
study the bulk conductance and also partition the system 
into three regions: the right (R) and left (L) leads - two 
semi-infinite gold mono-atomic chains - and a central (C) 
region, constituted by a single gold atom. This has the 
purpose of clarifying the role of both intra-conductor and 
conductor-lead correlations. We calculate the retarded 
Green's function in the space spanned by the MLWF set 
by inverting the Dyson equation, i.e. 



G''{uj) = [uj - H° -Y.^'iuj)]- 



(2) 



where if" is the KS hamiltonian once the exchange- 
correlation contribution is subtracted, i/'^ = Hks — Vxc. 
For the tri-partitioned geometry, the total retarded self- 
energy. 



(3) 



is the sum of the correlation GW and lead self-energies. 
The conductance is finally calculated using the following 
formula: 



G(e) 



2e^ 



tr[G"rHGTL(ri + rH)-ir[, 



(4) 



first given by Meir and Wingreen [l9| and recently re- 
derived under more general conditions 11|. Here F = 
Tl +ri^-|-rGWj not to be confused with the vertex func- 
tion, is the total decay rate (the sum of the electron in- 
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FIG. 2: (color online) (a) and (b): DFT-LDA Kohn-Sham 
vs GW electronic structure: black dots: DFT-LDA on PW 
basis; lines: DFT-LDA on MLWF; squares, triangles and dia- 
monds refer to GW calculations from the PP model, the CD, 
and the AC methods respectively. Fermi energy is set to zero, 
(c) Spectral function and (d) real and imaginary part of the 
GW self-energy: dashed, thin and thick lines are from the 
PP model, the CD, and the AC methods respectively. The 
dotted curve is the straight line uj — S + (14c) whose inter- 
sections with the real part of the self-energy give the peaks of 
the spectral function. The is set to the Kohn-Sham Fermi 
energy. 

and out-scattering functions), due to both the presence 
of the R and L leads and the effect of the e-e interaction, 

r = i(E''-I]'^) = i(S> -I]<), (5) 

With respect to the Landauer formula, Eq. ^ presents 
a factor (Fl -|- ri^)~^r which reduces to one in the un- 
correlated case. Equation ([4]) represents an effective (in- 
cluding correlation effects) transmission probability of an 
electron injected at energy e through the conductor. 

In Fig.[2i;a,b) we compare the DFT-LDA (Kohn-Sham) 
and the GW (quasiparticle) electronic structure for the 
relaxed geometry. The dots represent the DFT-LDA lev- 
els from the PW calculation. We verified that the diago- 
nalization on the MLWF basis (solid lines) closely repro- 
duces the PW results. Squares, triangles, and diamonds 
represent the GW electronic structures calculated using 
the PP model, the CD, and the AC methods, respec- 
tively. Little difference among the GW methods is found. 
The GW corrections globally lower the d-like states wrt 
the Fermi level, and also reduce the s-like bandwidth 
(this effect is less evident in the stretched chain where 
the bandwidth shrinks about 3 eV already at the DFT 
level). 

In Fig. [21 we compare the real and imaginary parts 
of the self-energy (d) and the spectral function A = 
i{G^ — G") (c) for a point close to the Fermi level. We re- 
mark that the AC method appears to smooth the richer- 
in-structure CD spectrum. Although computationally 
cheaper, the AC method should be considered less accu- 
rate than CD, especially on the imaginary part. However, 
the frequency dependence as well as the shape and the 
position of the main structures (both QP and satellites 
peaks) are essentially caught by both methods. There- 
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FIG. 3: (color online) Conductance (top) and spectral func- 
tion (bottom) for the stretched atomic configuration. The 
Fermi level is set to zero. Thin dashed line: Landauer result 
using a DFT KS structure; thin solid line: Landauer using 
only a real part GW renormalization of the energies; thick 
dashed line: Meir-Wingreen result using GW real part renor- 
malization in the leads and a full (hermitean-|-anti-hermitean) 
dynamical GW self-energy in the conductor; thick solid line: 
GW bulk conductance with full dynamical self-energy. 



fore, we use the AC approach in the following. 

In Fig. [3] we show the conductance and the spectral 
function of the gold chain for the stretched geometry, 
obtained using different methods: The thin dashed line 
is calculated using the Landauer formula and the DFT 
KS electronic structure. The conductance at the Fermi 
level for the stretched chain appears to be one (in units 
of 2e^//i) and it is of s-like character. This is true also for 
the relaxed structure (not shown), although in that case 
the Fermi level is at the limit of the onset of the d-like 
states. The thin solid line is obtained from the Landauer 
formula evaluated using the GW real-part-only QP en- 
ergies. GW corrections are considered both in the con- 
ductor and in the leads. Otherwise, a ficticious contact 
resistance, unphysical for a homogeneous system, would 
appear. At this first level, the net effect is a renormaliza- 
tion of KS into QP energies, the true energies to intro- 
duce and remove an electron from the system. Therefore, 
the GWA affects the conductance profile by modifying 
the position of the conductance steps, especially in the 
d-like region. In the relaxed geometry, the GWA also 
narrows the s-like conductance channel. 

The thick dashed line in Fig. [3] represents the result ob- 
tained in the tri-partitioned geometry by using the Meir- 
Wingreen formula and introducing a full non-hermitean 
and dynamical GW self-energy in the conductor. Static 
real-part-only QP energies are included in the leads. This 
introduces loss-of-coherence only in the conductor while 
leaving the leads ballistic. At the same time it limits 
the introduction of ficticious contact resistances, i.e. the 
QP levels are aligned in the leads and the conductor. 



4 



The difference of this curve wrt the thin sohd hne gen- 
uinely represents the effect of e-e scattering mechanisms 
in the conductor, causing diffusion, loss of coherence and 
appearance of resistance. With respect to Landauer ap- 
proaches, the spectral function now appears as a collec- 
tion of broadened QP peaks, whose finite width is di- 
rectly associated to the inverse of the electronic lifetime 
of the QP state. The spectral weight, which is spread 
out, results in a lowering and a spill-out of the conduc- 
tance step-like profile. This effect is directly related to 
the imaginary part of the QP energies, and can be seen 
to increase with e — Ep although not with a quadratical 
scaling Fermi-liquid behaviour, as it is normally observed 
in GW results for 3D systems. Finally, we calculate the 
fully correlated bulk GW conductance by taking into ac- 
count a non-hcrmitian and dynamical GW self-energy ev- 
erywhere in the system, conductor and leads (thick solid 
line). With respect to the previous case, even residual 
contact resistances (due to the fact that the conductor 
and leads spectral peaks were differently shaped, with fi- 
nite and infinitesimal widths respectively) are completely 
removed, and the conductance increases almost overall. 
Only around —3 eV we see a slight drop, which is due 
to the specific (Fl -|- Ffl;)~iF factor in Eq. (g]). More- 
over, new structures appear in the conductance at the 
lowest energies. By inspecting the spectral function, we 
can attribute them to the presence of satellites of elec- 
tronic origin, i.e. plasmons or shake-ups, of the main QP 
peaks. Since the e-e interaction is an elastic scattering 
mechanism, these satellites are necessary to balance the 
losses which occur at energies close to the Fermi level, and 
are therefore important for transport. The e-e scattering 
acts in a way to redistribute the conductance channels to 
different energies, rather than globally destroy conduc- 
tance as in the e-ph scattering, where momentum and 
current is lost to ionic degrees of freedom. 

Taking the GW bulk result, we have integrated the 
conductance curve such to obtain the voltage character- 
istics of the correlated system. We compare with the 
experimental results of Ref. [2^ and the e-ph result of 
Ref. [§], calculated at exactly the same stretched 5.35 
Bohr interatomic distance. Like in that work, we as- 
sume that an equilibrium picture can still be appropriate 
to describe the small voltage range of ± 30 mV. The 
GW result is shown in Fig. 2] (dots). The results from 
Ref. attribute the step in the conductance, occurring 
at ~ 15 mV, to the onset of phononic processes. Instead, 
the continuous drop observed in our electronic correlated 
conductance, occurring in the first 15 mV, compares fa- 
vorably with the drop observed experimentally [20|: e- 
e scattering mechanisms seem hence responsible for the 
conductance drops at very low bias. While the quantita- 
tive agreement with the experiment on the conductance 
value may be somewhat fortuitous [27| , the trend in this 
drop is a direct consequence of the increase in the GW 
imaginary part of QP energies. 
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FIG. 4: (color online) Diflterential conductance vs applied 
bias. Thin solid line: DFT Landauer result; dots: present 
NEGF GW bulk result for the 5.35 Bohr interatomic distance; 
dashed and dotted line: e-ph theory of Ref. [!|] correspond- 
ing to 4 atoms, same interatomic distance and for the damped 
and undamped limits; thick solid lines: experimental result of 
Ref. [2^ corresponding to 2 and 7 atoms and different chain 
strains. 



In conclusion, we have calculated the conductance of 
a realistic gold chain system by taking into account e-e 
correlation effects within the GW approximation. With 
respect to Landauer DFT results, the conductance pro- 
file is considerably modified. Already at the level of an 
equilibrium non-self-consistent GW, the trend of the dif- 
ferential conductance appears to compare favorably with 
the trend experimentally observed for this system. 
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